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We review recent progress in calculating the twist— 2 parton distribution functions 
(PDF's) at a low normalization point in the large-iV c limit, where the nucleon can be 
CN \ described as a soliton of the effective chiral lagrangian. This field-theoretic approach 

preserves all general requirements on the PDF's (positivity, normalization etc.). In 
. particular, it allows to calculate the polarized and unpolarized antiquark distributions 

! at the low normalization point. 

In this talk we review recent attempts to calculate the leading-twist parton distributions 
of the nucleon at a low normalization point in the large- N c limit [|], |2|, H|. The results 
presented here have been obtained in collaboration with D.I. Diakonov, V.Yu. Petrov, P.V. 
Pobylitsa, and M.V. Polyakov. 
^ Introduction. The parton distribution functions (PDF's) of the nucleon are an essen- 

tial ingredient in the description of deep-inelastic scattering and a variety of other hard 
processes. The scale dependence of the PDF's in the asymptotic region is determined by 
the renormalization group equation of perturbative QCD. A complete description of the 
experiments requires non-perturbative information in the form of the values of the parton 
distributions at some initial normalization point. Several sets of input distributions have 
been determined by parametrizing data from a number of different experiments at large q 2 
[|], [51 ^ 0| . All modern parametrizations involve both antiquarks and gluons at the low scale. 

When attempting to compute the parton distributions at a low normalization point from 
first principles one must deal with non-perturbative effects such as the dynamical breaking 
of chiral symmetry and confinement. While a satisfactory theory of confinement is still 
lacking, the dynamical breaking of chiral symmetry in QCD is well understood theoretically: 
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the resulting effective theory at low energies is known, and a dynamical picture of chiral 
symmetry breaking is provided by the instanton vacuumF] []T0|, |ll[] , which has recently received 



direct support from lattice simulations [|12}| . A description of the nucleon based on the 
dynamical breaking of chiral symmetry is possible in the large-iV c limit, where QCD becomes 
equivalent to an effective theory of mesons, and the nucleon emerges as a soliton of the pion 
field |TB|, |TJ], [T5[ . This picture of the nucleon is known to give a good account of hadronic 



observables such as the nucleon mass, magnetic moments, form factors etc. |16[. It is 
thus natural to describe also the parton distribution at a low normalization point in this 
approach. An important advantage of this description of the nucleon is its field-theoretic 
character, which is essential for the parton distributions to satisfy general requirements such 
as positivity, proper normalization etc. 

Effective chiral theory. In the large- N c limit, the consequences of the dynamical breaking 
of chiral symmetry can be encoded in an effective lagrangian for the pion field, valid at low 
energies. It can be expressed as a functional integral over quark fields (SU(2) flavor) in the 
background pion field, tt(x), |TT|, [T7|, |TS 



exp(iS cS [n(x)}) = J DipDtfj exp i J d^x^iitf - MIP 5 )^ , (1) 
W 5 (x) = exp [ivr a (x)r a 7 5 ] . (2) 

Here, M is the dynamical mass (generally speaking, it is momentum-dependent), which is 
due to the spontaneous breaking of chiral symmetry. Eq. ([]]) describes the minimal chirally in- 
variant interaction of quarks with Goldstone bosons. In the long-wavelength limit, expanding 
in derivatives of the pion field, the effective action Eq.(Q) reproduces the Gasser-Leutwyler 
lagrangian with correct coefficients, including the Wess-Zumino term. The effective theory 
defined by Eq. (|TJ) is valid for momenta up to an UV cutoff, which is the scale at which the 
dynamical quark mass drops to zero. For simplicity we shall take in the discussion here 
the quark mass to be momentum-independent and assume divergent quantities to be made 
finite by applying some regularization scheme; below we shall see why, and under which 
conditions, this approximation is justified when computing (anti— ) quark DF's. 

The effective action Eq.(|l]) has been derived from the instanton vacuum, which provides 
a natural mechanism of dynamical chiral symmetry breaking and enables one to express the 
parameters intrinsic in Eq.(|T]) — the dynamical mass, M, and the ultraviolet cutoff — in 
terms of the QCD scale parameter, Aq C d [PH - The cutoff is given by the average instanton 



size, 

p- 1 ~ 600 MeV, (3) 

which is the order of the scale at which the coupling constant is fixed in the instanton vacuum. 
Since the effective theory is not valid beyond p _1 , when applying it to compute PDF's we 
are certainly dealing with distributions normalized at a scale not higher than 600 MeV. 

An important point to emphasize is that, when working with the effective action, Eq. ([!]), 
it is implied that the ratio of the dynamical quark mass to the ultraviolet cutoff is paramet- 
rically small, Mp < 1. In the instanton vacuum this is guaranteed by the fact that M 2 is 

1 For a recent review of the instanton vacuum sec refs.tfsl M. 
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proportional to the instanton density, so that {Mp) 2 is proportional to the packing fraction 
of the instanton medium, 



{Mpf oc (^j , (4) 

which is a small parameter. However, this requirement is really of a more general nature - 
otherwise degrees of freedom other than those included in Eq. ([!]) would play an essential role 
in the effective dynamics, and the concept of an effective action would not be very meaningful. 
One may interpret the inverse ultraviolet cutoff, p, as the size of the "constituent" quark 
described by the effective theory. The smallness of Mp implies that the constituent quark is 
nearly pointlike — its size is small compared to its typical wavelength in the nucleon, which 
is of order] 2 ] M _1 . This fact is essential for the interpretation of the distribution functions 
computed with the effective chiral theory. For Mp —>■ one may identify them with the QCD 
quark and antiquark distributions. When trying to be more accurate, keeping terms of order 
(Mp) 2 , one begins to resolve the structure of the "constituent" quark, and the distribution 
functions computed in the effective theory should be regarded as distributions of composite 
objects which have a substructure in terms of QCD partons. 

To describe the gluon distribution one needs a microscopic picture of the non-perturbative 
gluon degrees of freedom. Such a picture is provided by the instanton vacuum, which on one 
hand allows to derive the effective chiral action, Eq.([TJ), by explicit integration over gluon 
degrees of freedom, on the other hand makes it possible to directly evaluate matrix elements 



of gluon operators, using the effective operator method developed in ref.|19[. One finds that 
the twist-2 gluon distribution is of order (p/ R) A , which is consistent with the above statement 
that the structure of the constituent quark starts to appear at level (Mp) 2 . To compute the 
twist-2 gluon distribution one needs to construct both the effective action and the effective 
operators at order (Mp) 2 . In the computation of twist-2 quark and antiquark DF's reported 
here we work in the limit Mp — > 0, where the constituent quark is pointlike and the twist-2 
gluon distribution is zero. For a more detailed discussion of the gluon distribution, including 
the role of the gauge field in the twist-2 quark operators, see refs. ]2T], [22| .[] 

The nucleon as a chiral soliton. In the effective theory defined by Eq.(|l|) the nucleon is 
in the large- N c limit described by a static classical pion field ("soliton") ||15|| . In the nucleon 
rest frame it is of "hedgehog" form, 

f/ c 75 (x) = exp [m a r a P(r) 75 ] , r = |x|, n = ~' ( 5 ) 

where P(r) is called the profile function, satisfying -P(O) = — ir, P{r) — > for r — > oo, which 
is determined by minimizing the energy of the static pion field. Quarks are described by 
single-particle wave functions, which are the solutions of the Dirac equation in the external 



2 The dominant contribution to the twist-2 quark and antiquark DF's comes from quarks with momenta 
of order M <C see below. 

3 Calculations of the twist-2 gluon distribution in the instanton vacuum using a different approach have 
recently been performed by Kochelev pfl]. However, the contributions taken into account there do not 
represent the full answer to order {p/R) , see refs. pi], |22^] for a discussion. 
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pion field, 



H$ n = E n $ n , H = - n V4 + M 7 ^. (6) 

The spectrum of the one-particle Hamiltonian, H, contains a discrete bound-state level, 
which must be occupied by N c quarks to have a state of unit baryon number, as well as the 
positive and negative energy Dirac continuum, distorted by the presence of the pion field. 
The nucleon mass is given by the minimum of the bound-state energy and the aggregate 
energy of the negative Dirac continuum, the energy of the free Dirac continuum subtracted 
[|Tj|. Nucleon states of definite spin-isospin and 3- momentum are obtained by quantizing 
the rotational and translational zero modes of the soliton, i.e., applying a flavor rotation 
and a shift of the center to the soliton field, Eq.flU), and integrating over the corresponding 



collective coordinates with appropriate wave functions [15, |i6fl. 



Quark- and antiquark distributions from the effective chiral theory. In order to compute 
the twist-2 quark and antiquark distribution functions in the effective chiral theory one may 
start either from their "parton model" definition as numbers of particles carrying a given 
fraction of the nucleon momentum in the infinite-momentum frame, or from the "field- 
theoretic" definition as forward matrix elements of certain light-ray operators in the nucleon. 
Both ways lead to identical expressions for the quark DF's in the chiral soliton model. This 
remarkable fact should be attributed to the field-theoretic nature of this description of the 
nucleon, and to the fact that the main hypothesis of the Feynman parton model — that 
transverse momenta do not grow with q 2 - - is satisfied within this model. In the first 



formulation one considers a nucleon with large 3-momentum, 



and defines the quark and antiquark DF's as 



M N v 

P n = H u -> 1, (7) 




d z h ( 1.3 \ f </(kK/(k) 

;2n5 { x - ^) (N v \ \ } |iV v >, (8) 



(2tt) 3 V P, 



N 



where \N V ) denotes the state with the fast-moving nucleon and afj,a it f and bfj,bij are 
the creation/annihilation operators for quarks (antiquarks), characterized by a polarization 
quantum number i and flavor /. For large N c the nucleon consists of 0(N C ) quarks and 
antiquarks, so we can hope to describe the distribution functions for values of x of order 

In the large N c limit the nucleon matrix element Eq. (|8]) reduces to a sum of matrix elements 
between quark single-particle states in the mean pion field (saddle point approximation), 
which is given in this case by the hedgehog field, Eq.(||), boosted to velocity v, see || for 
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details. Taking the limit v — > 1 one arrives at the following expressions for the twist-2 quark 
and antiquark distributions of the nucleon in the large- N c limit: 



D 



i,f 



X 



A,/(x) 



N r M 



c^N 



E 

occ. 

E 



non— occ. 



where, for example, 



/ 72^^ (k)(1 + T°7 3 )7°IV(fc 3 + E n T xM N )$ n>f (k), 



(2 



7 



,1±75 



(10) 



(11) 



for quarks polarized along or against the direction of the nucleon velocity. Here the sums run 
over all occupied or non-occupied quark levels in the hedgehog pion field, that is, the bound- 
state level and the negative Dirac continuum or the positive Dirac continuum, respectively. 
It is understood in Eq.([T0|) that one subtacts the corresponding sum over eigenstates of the 
free Hamiltonian, where ir(x) = 0. 

Alternatively, one can start from the QCD definitions of the twist-2 quark distribution 
functions as matrix elements of certain light-ray operator, which act as generating functions 
of the series of local twist-2 operators [E4L Bf| ; in the unpolarized case 



-^unpol,/ (*^) 
-^unpol,/ (•^) 



- j dz-e^*-(p\MQ)i + M*)\P) 

— CO 



2+=0, x± =0 



{x — > —x} 



and in the longitudinally polarized case 



(12) 



1 / . 



Dpo\j(x) = — I dz e 



47T 



D po lj{x) = {X -> -x} . 



(P,5|^ / (0) 7 + 75^/(^)|P,^) 



Z+=0, 2 I =0 



(13) 



These matrix elements can be evaluated in the nucleon rest frame, expanding the quark fields 
in Eqs.(|T2"|, fTS] ) in the basis of single-particle wave functions, Eq.(^); the result is identical 
to Eq.flnj). 

To obtain the unpolarized or polarized (anti-) quark distributions for a nucleon state of 
definite spin and isospin, one has to apply a flavor rotation to the quark wave functions in the 
basic expression Eq . (|T0|) , $ n j — > Rff>(t)§ n j', and integrate over slow rotations, projecting 
on a nucleon state with given spin and isospin ffl. At this point characteristic differences 
between the different isospin combinations of the unpolarized and polarized distributions 
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appear. The angular velocity of the soliton, F$(dR/dt), is of order 1/N C (the soliton moment 
of inertia is 0(N C )), and in the leading order of the l/iV c -expansion one may neglect it. 
The distribution functions which are non-zero at this level are the isosinglet unpolarized 
and isovector polarized — they are the leading ones in the l/iV c -expansion. The isovector 
unpolarized and isosinglet polarized distributions, on the other hand, are non-zero only after 
expanding to first order in the angular velocity; they are thus suppressed relative to the 
former by a factor 1/N C . Thus the l/iV c -expansion leads to the following classification of 
quark/antiquark DF's in "large" and "small" ones: 





unpolarized 


polarized (longitud.) 


polarized (transv.) 


"large" 


u + d 


Au - Ad 


h\ - h\ 


"small" 


u — d 


Au + Ad 


h\ + hf 



More precisely, bearing in mind that x 
the large N c limit behave as 



£)large / 



N?f(N c x) 



1/N C , we can say that the (anti-) quark DF's in 
D smal \x) ~ N c f(N c x), (14) 



where f(y) is a stable function in the large iV c - limit which depends on the particular distribu- 
tion considered. One may easily convince oneself that this is consistent with the well-known 
^-dependence of the lowest moments of the DF's: the first moments of u + d (baryon num- 
ber) and of Au — Ad (isovector axial coupling, g^) are 0(N C ), while those of u — d (isospin) 
and of Au + Ad (isosinglet axial coupling, g^') are 0(N®). 

The "large" DF's are given by simple sums over quark levels; the isosinglet unpolarized 
DF by (the sum over flavor indices is understood on the R.H.S) 



— — d>t (k)(l + 7 V)5(fc 3 + E n + xM N )$ n (k), (15) 
[u + d)(x) = N C M N V {x^-x}. (16) 



An alternative representation of the DF is obtained by applying in Eq. (|l2|) the anticommuta- 
tion relation of quark fields before taking the light-cone limit; this leads to a representation of 
the quark DF's as sums over non-occupied, and of the antiquark DF as a sum over occupied 
states: 



[u + d](x) 



-N C M N J2 {x^-x} 



(17) 



where the braces again denote the corresponding expressions appearing in Eq.(|15]). One may 
regard Eq.(|l5|) together with Eq.([L7]) as one universal function, describing the quark DF at 
x > and minus the antiquark distribution at x < 0. The isovector-polarized DF is given 
by 4 

1 r d 3 k 

[Au-Ad)(x) = - -(2T 3 )N C M N Y: /-— I $J i (k)(l + 7 7 3 )75^ 3 + J E„ + xM^)$ n (k), 

& occ. J {^ 7r ) 



For the corresponding formulas in the case of transverse polarization, see ref. | 
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(18) 

[Au-Ad](x) = ~(2T 3 )N C M N £ {x - -x} , (19) 

non— occ. 



or, alternatively, 



[Au-Ad}(x) = -l-(2T 3 )N c M N J2{x^-x}, (20) 

" occ. 



where 2T 3 = ±1 for proton and neutron, respectively. The "small" DF's in the large-iV c 
limit are given by double sums over quark levels divided by the moment of inertia of the 
soliton; see ref. |J for details. 

From the explicit representations Eqs.(|T5|-p(]|) it can easily be seen that the large- N c 
quark and antiquark DF's satisfy the following properties: 

Positivity: The singlet unpolarized quark and antiquark DF's, Eqs.([T^, [TB|), are explicitly 
positive (the Dirac matrices in these expressions are a projector and thus positive definite). 
Sum rules: Integrating Eq. (|i~5|) minus the representation Eq. ([l7|) over x one finds 



J dx[u + d — u — d]{x) = N c (number of occupied levels in the soliton 

— number of occupied levels in the vacuum) , (21) 

which coincides with the baryon number of the nucleon in the large- N c limit. One can also 
show that the momentum sum rule is satisfied by virtue of the fact that the soliton field is 
a stationary point of the static energy. Similarly, integrating Am — Ad + Am — Ad, Eq. (IT$) 



(3) 

minus Eq.(20), one obtains the expression for the isovector axial coupling, g A , in the chiral 



soliton model |T5], [16[ , i.e., the Bjorken sum rule is satisfied within this approach. 

For a proof of sum rules for the "small" DF's (e.g. the isospin sum rule for u — d), and a 

discussion of the Gottfried sum we refer to ref.]l[. 

It is essential that both the contributions of the bound-state level as well as the Dirac 
continuum are taken into account in the representation of the large-iV c quark and antiquark 



DF's, Eqs.(|15]-|2Cj). Restricting oneself to the bound-state level only one would, for example, 
violate either the the baryon number sum rule or the positivity of the singlet unpolarized 
antiquark distribution — since the level makes a negative contribution to Eq . (|17|) , depending 



on whether one writes the antiquark distribution in the form Eq.(]TB|) or Eq.([T7]). Adding the 



contribution of the Dirac continuum Eq.(^) and Eq.(|P7D give identical results, the antiquark 
distribution is positive, and the baryon number unity. 

In the discussion of the large- iV c DF's so far we have not explicitly taken into account the 
UV cutoff, which is an essential ingredient in the effective chiral theory, Eq.(|I|). In fact, the 



expressions for the distribution functions as sums over quark levels, Eqs.(|l5| - ^0|) contain 
UV divergences which are made finite by the UV cutoff. When choosing a regularization 
method to implement this UV cutoff it is of utmost importance that the regularization do 
not violate any of the fundamental properties of the DF's. A crucial requirement is that 
the regularization should preserve the completeness of the set of quark single-particle wave 
functions in the soliton field, Eq.(O). An incomplete set of basis functions would mean a 
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violation of the local anticommutation relation of the quark fields, that is, a violation of 
causality. For an extensive discussion of the many facets of this we refer to ref.@. Let 
us mention here only that the equivalence of the representations of the (anti-) quark DF's 
as sums over non-occupied or occupied states, cf. Eqs. (|16|, |TTD , which is instrumental in 



ensuring both positivity and proper normalization of the isosinglet DF, depends crucially on 
the completeness of the basis. 

A regularization by subtraction, for example, a Pauli-Villars cutoff, meets the above 
requirement, while a regularization by cutoff (for example, an energy cutoff or the popular 
proper-time regularization of the fermion determinant) is equivalent to working with an 
incomplete set of states and thus unacceptable. The Pauli-Villars cutoff is implemented by 
subtracting from the sums Eqs.(^ - |20|) a multiple of the corresponding sums computed 
with a regulator mass, M PV ~ 



D(x) PV = D(x) 



M 



2 



M 2 

M IV1 PV 



D(x) 



(22) 

M PV 



It can be shown that such regularization preserves all qualitative properties of the DF's 
(positivity, sum rules) for which the completeness of states is an essential condition, thus 
justifying our above discussion of "naive" unregularized expressions 0. (Note that only 
UV divergent quantities such as the total isosinglet distribution u + d + u + d should be 
regularized; UV finite quantities such as the valence distribution u + d — u — d should be left 
unregularized.) In short, the UV cutoff in the effective theory Eq. (jlj) can be implemented 
in a way that the field-theoretic character of the description of quark and antiquark DF's is 
preserved. 

It is important to note that, with a regularization preserving completeness, the actual 
dependence of the twist-2 DF's on the cutoff of the effective theory is very weak: the DF's 



Eqs.(|l5| - 120|) are only logarithmically divergent for fixed x, and so are all their moments, 
see ref. [||. The dominant contributions to the DF's thus come from degrees of freedom 
with momenta of order M <C This fact is crucial for the consistency of the effective 
theory approach: the quark and antiquark DF's are determined by genuine dynamics of 
the effective theory — that is, by degrees of freedom with momenta much smaller than 
the cutoff — not by the details of the UV regularization, which, at any rate, is known 
only schematically. In particular, this justifies working with a suitable generic regularization 
(Pauli-Villars subtraction) rather than with the momentum-dependent quark mass obtained 
from the instanton vacuum. 

One exception to the above is the region of values of x < (Mp) 2 /N c , where the dominant 
contribution to the (anti-) quark DF's comes from states with large (positive or negative) 
energies^. In this parametrically small region of x the momentum dependence of the quark 
mass is essential and, in the case of the singlet unpolarized distribution, Eqs. (|T5|, |TTD , leads 
to a vanishing of the distribution at x = 0, as has recently been shown in connection with 



5 We remind that we are discussing the (anti-) quark DF's at a low normalization point, whose small-a; 
behavior has no direct relation to that of the DF's at experimental scales, which is determined essentially 
by perturbative evolution. 
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xD(x) 




Figure 1: Solid lines: The isosinglet unpolarized valence quark and antiquark distributions 
computed in the large-iV c limit Jl|, §]. Dashed lines: NLO-parametrization of GRV ||. 



an investigation of the off-forward PD's in the large- iV c approach [27]. We note also that 



in contrast to the twist— 2 DF's, matrix elements of higher-twist operators generally receive 
contributions from momenta of order of the cutoff, p _1 , when computed in the effective 
theory. Such matrix elements can be estimated within the effective theory derived from the 
instanton vacuum, taking into account the full momentum dependence of the constituent 
quark mass [EI], EES. 

Results for quark and antiquark distributions. To compute the DF's one needs to eval- 



uate the sums over quark levels, Eqs.(|T5] - E_Up, taking into account the contributions of 
the bound-state level as well as the polarized Dirac continuum. This can be done either 
by numerical diagonalization of the hamiltonian f2j, or by rewriting the sums over levels in 
terms of the quark Green function in the background pion field, which allows an approximate 
evaluation of the continuum contribution using the so-called interpolation formula . 

Results for the isosinglet-unpolarized and isovector-polarized distributions are shown 
in Fig.l and 2, where we compare with the parametrizations of Gliick, Reya et al. || 0. 
In the calculations reported here we have chosen a simple UV regularization by way of 
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Figure 2: Solid lines: The isovector polarized (longitud.) total distribution (quarks plus 
antiquarks) and antiquark distribution computed in the large- N c limit |], |2|. Dashed lines: 
NLO-parametrization of ref.0. The polarized antiquark distribution is assumed zero in the 
fit. 

Eq.fl22|) and used the standard parameters for the effective chiral theory (see refs.|], @ for 
details); we have not attempted to make a best fit. Fig.l shows the isosinglet unpolarized 
valence and antiquark DF. We note that the antiquark distribution obtained in our approach 
is of the "valence-like" (i.e., non-singular at small x) form assumed which was assumed 
in the fits of ref.0. Fig. 2 shows the isovector polarized total distribution (quarks plus 
antiquarks) as well as the antiquark distribution, which was set to zero in the fit of ref . |F] . 
That the calculated total distribution is systematically smaller than the fit is due to the 
fact that the isovector axial coupling, is underestimated by the large-iV c limit; the 



1/N C — corrections in this channel are large and have been computed [16]. In Figs.l and 2 we 
compare the calculated DF's with the parametrizations at face value, without taking into 
account evolution; the normalization point of the calculated DF's may actually be lower 
than that of the parametrizations. 

We stress that we are computing the twist-2 parton distributions at a low normalization 
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point, not the structure functions (cross sections) at low q 2 , which are, in principle, directly 
measurable. The latter are affected by higher-twist (power) corrections, which become large 
at low q 2 . Our calculated distributions should be used as input for perturbative evolution, 
starting with a scale below the cutoff, p _1 ~ 600 MeV. A direct comparison with the data 
can be performed only after perturbative evolution to sufficiently large q 2 . 

Results of a calculation of the "small" DF's in the large- N c limit, the isovector unpolar- 
ized and isosinglet polarized, have been reported by Wakamatsu and Kubota [^6] . 



Recently also the off-forward quark DF's have been computed in the chiral soliton model 
27fl . These generalized DF's are defined as nonforward matrix elements of light-ray opera- 



tors of the type of Eq. (non-zero momentum transfer) . It was found that the large-iV c 
approach predicts a strong dependence of the isosinglet off-forward distribution on the lon- 
gitudinal momentum transfer. 

Summary. To summarize, we have shown that the large-iV c picture of the nucleon as a 
soliton of the effective chiral theory is able to explain the essential features of the quark and 
antiquark DF's at a low normalization point. As we have shown, the success of this approach 
is essentially due to two circumstances: 

• the parametric smallness of the ratio of the dynamical quark mass to the cutoff of the 
effective theory, Mp (related to the small packing fraction of the instanton medium), 
which allows to identify the twist— 2 QCD (anti-) quark DF's with the distributions 
computed in the effective theory; 

• the field-theoretic character of this description, which ensures that the DF's satisfy all 
general requirements such as positivity, normalization etc. 

An important direction for future work is the refinement of this approach to a level that 
allows one to take into account effects of higher order in (Mp) 2 , i.e., to resolve the structure 
of the "constituent" quark. This can be done in the framework of the instanton vacuum, 
which allows to derive the effective chiral action by explicit integration over non-perturbative 
gluon degrees of freedom. In particular, at order (Mp) 2 the gluon distribution will start to 
appear. Also, this will allow to determine the normalization point of the calculated DF's 
more accurately. 
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